Discrete breathers in BEC with two- and three-body interactions in optical lattice 
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We investigate the properties of discrete breathers in a Bose-Einstein condensate with three-body 
interactions in optical lattice. In the tight-binding approximation the Gross-Pitaevskii equation with 
periodic potential for the condensate wavefunction is reduced to the cubic-quintic discrete nonlinear 
Schrodinger equation. We analyze the regions of modulational instabilities of nonlinear plane waves. 
This result is important to obtain the conditions for generation of discrete solitons(breathers) in 
optical lattice. Also using the Page approach, we derive the conditions for the existence and stability 
for bright discrete breather solutions. It is shown that the quintic nonlinearity brings qualitatively 
new conditions for stability of strongly localized modes. The numerical simulations conform with 
the analytical predictions. 
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I. INTRODUCTION 

Discrete breathers are important eigenmodes of the 
nonlinear discrete and periodic systems and their ex- 
istence has been predicted for such systems as optical 
waveguide arrays,array of Josephson junctions, magnetic 
chains etc [H El El- Recently discrete breathers in ar- 
ray of Bose-Einstein condensates has attracted a great 
deal of attention 0, IE B H H H EJ. The localized 
atomic matter wave in BEC in the optical lattice (a gap 
soliton) has been observed recently in ^l|. The theo- 
retical approaches mainly based on the analysis of the 
Gross-Pitaevskii equation with periodic potential and cu- 
bic mean field nonlinearity. It corresponds to taking into 
account elastic two-body interactions and the potential 
of optical lattice. The progress with Bose-Einstein con- 
densates on the surface of atomic chips and in atomic 
waveguides involves strong compression of BEC and so 
to essential increasing the density of BECs. Then it pose 
the problem of the taking into account three-body inter- 
actions effects. This interaction represent interest also 
for the understanding of the fundamental limits for the 
functioning of devices using BEC |l2j. The existence 
of three-body interactions can play important role in the 
condensate stability Recently 3-body interaction 

arising due to the Efimov resonance has been observed 
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in an ultra-cold gas of cesium atoms |15| . 

Thus represent interest to investigate properties of dis- 
crete breathers in BEC in optical lattices when two- and 
three-body interactions are taken into account. One of 
the interesting properties of this system is the existence 
of the gap-Townes soliton pl|. In the tight-binding ap- 
proximation(see a section below) we show that the GP 
equation with periodic potential can be reduced to the 
cubic-quintic discrete nonlinear Schrodinger(CQDNLS) 
equation. This regime occurs for the deep optical lattice, 
which is realized in typical experiments in the BEC. 

For the generation of discrete breathers the modula- 
tional instability of nonlinear plane wave represent the 
possible mechanism. Thus in this paper in the section 3 
we discuss the modulational instability in the CQDNLS 
equation, the stability region of the wave in the lattice. 
In section 4 we discuss the localized modes of the different 
symmetries and the conditions of their stability. 



II. THE MODEL 

Let consider BEC with two and three body interac- 
tions in optical lattice and in highly anisotropic trap cor- 
responding to quasi ID geometry. The Gross-Pitaevskii 
equation is 



m t = -— V 2 ^ + y (loIp 2 + uj 2 x x 2 ) + V opt (p, x) - ffilVlV - ffaMV, 



(1) 



r 



where g\, gi are coefficients proportional to two- and is expected that the three-body interaction is attractive 
three-body elastic interactions. For rubidium atoms it and \gs\ ~ 10 -26 — 10~ 27 cm 6 /s 12] . Here uij_ uj x , 



2 



V op t = Vq sin 2 (kx)R(p). We look for the solution of the 
form 

il>(r,t) = <fo(p)<Kx,t), (2) 
where <po is the solution of the radial linear equation 

h 2 2 m 2 2 1 
— V p o + 00 = (3) 

The solution for the ground state is 



Multiplying both sides of the GP equation by 4>o and 
integrating over transverse variable p we obtain the quasi 
ID GP equation 



? eXp( -2^ } - 



itujn = -^—<f)xx + Ctt^Ix 2 + V a sm 2 (fc:r))</> 
2m 2 



,'/2 



37r 2 a 



(4) 



Introducing the dimensionless variables 



t = tv, x = kx, 4> = u 



2q s uj± 



-,v = 



Er 



where Er = h 2 k 2 /2m is the recoil energy, we obtain the 
equation 



' \\u\ 2 u + \\u\ 4 u- V cos(2x)u = 0. (5) 



iu t + -u xx + 



where A = ±1 and 
A 



!J2 



Z~K 2 fwa\ \ 2a a w± 
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In the case of deep optical lattice (U 3> Er i.e. c>l) 
the tight-binding approximation can be used 0, El IE 0] • 
If the potential U(x) periodic, i.e. U(x + I) — U(x), the 
eigenvalue problem is: 

—T^T +U(x)(t> a ^ q = E a ,q<t>a,q, (6) 

where for the optical lattice potential U(x) = Vq cos(2x). 
Here <f> a ,q is the Bloch function, a is the index labelling 
the energy bands E a (q) . Periodicity of the lattice admits 
the expansion of the energy 



E, 



'a, 9+90 



— E aA — E a - q , qo — 



2ir 

T' 



Thus we can expand the energy in the Fourier series 

I rio/2 



E a . q - 2J 



lqnl '•' = — / E q ^e- lqnl dq. 



2tt J_ 



9o/2 



(7) 



For deep lattice case it is convenient to use the Wannier 
functions 0: 



I f qa/2 



I TU,9 

27r -/-90/2 



a , q (x)e-^ nl dq, (8) 



and 



2tt 



w n<a {xy« nl . 



(9) 



The main property of Wannier functions that they are 
strongly localized in the bottoms of the potential and 
so are suitable for the description of strongly localized 
modes in the periodic potential. They form an orthonor- 
mal and a complet set 



^n^ce^n' ,oc' \Xjdx — 3n,n' &a.a f : 

> 

*yi w n,a( x ) w n,a(z') = 8{x - x'). (10) 



Let look for the solution of CQ GP equation of the form 
of expansion in Wannier functions 



u(x, t) 



(11) 



Substituting this expansion into the GP equation JSJ) we 
obtain the system of equation for coefficients c n , a 
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i,a + A 



E 



E c ^ c 3^<= + A E E 4c* 2 c 3 c 4 c 5 Wt, 



(12) 



5a 5n 



r 



where 



nni n2ri3 



For deep optical lattice Vq > 5Er are important terms 
with n = n\ — — W53 i.e. we are restricted by the 
band a. Then the system is reduced for equation 



(C n _i lQ + G„_i, a ) + A^|C„| 2 C„ + AW 6 ^|C„| 4 C„. 



(13) 



r 



By the transformation v n ■ 
can be transformed into 

iv n ,t + n{v n +i + V n -i) + 



Cn,a\/Wf^Cn this equation 



A|v„| 



0. (14) 



Here k = wf jCt , 7 = \W$£ / '(2(Wf™) 2 ) . It is the cubic- 
quintic discrete nonlinear Schrodmgcr(CQDNLS) equa- 
tion. As showed the analysis performed in Q , the region 
of parameters exists, when the tight-binding approxima- 
tion describes the GP equation with periodic potential 
with a high accuracy. It is achieved for the amplitudes of 
the periodic potential |Vq| > 5-Er, and Er is the recoil 
energy. That corresponds to the experimentally realized 
values. When 7 = the problem is reduced to considered 
previously by 0, 0, ■ In the recent work 01 the case 
when A > 0, 7 < 0, motivated by the analysis of localized 
states in array of waveguides with saturable nonlinear- 
ity, has been considered. In the BEC case the parameter 
7 can has any sign. The conserved quantities arc: The 
norm(number of atoms) N is: 



^ = Eki 2 ' 

n 

and the Hamiltonian H is: 

E 

n— — c 

The equation of motion is: 



(15) 



n{v* n v n+ i + v n v* +1 ) + -K| 4 + ^K| 6 



5H 

Su* 



(16) 



(17) 



The momentum is not conserved due to the discretness 
of the system. Due to the periodicity of the system the 
Hamiltonian is periodic with period n = 1. The differ- 
ent excitations induces different effective periodic poten- 
tial like the Peierls-Nabarro potential. The height of the 



barrier depends on the localization of the discrete soli- 
ton solution - grows with the decreasing of the soliton 
width. The analysis of the Hamiltonian can help for the 
calculation of the PN barrier between different discrete 
solitons and so help for the investigation of the stability 
of nonlinear modes 0] . 



III. MODULATIONAL INSTABILITY OF 
NONLINEAR PLANE WAVES 



The nonlinear discrete equation has the plane wave 
solution 



v n (t) = u e^ n+ut \ 
with the dispersion relation 

ui = 2k cos q + A^o 2 + 7^o 4 - 



(18) 



(19) 



For special importance here are the staggered q = tt 
and unstaggered q = solutions 0] where the dispersion 
relation can be rewritten as 



u = ±2k + \v 2 + 7^0 4 



(20) 



The linear stability of the solution (|18fl and (|19fl can be 

investigated by looking for a solution in the form 



V n (t) = (u Q +5v n {t)y {qn+ ^\ 



(21) 



where Sv n (t) is a small perturbation on the carrier 
wave, substituting (|21|l in (|14l) and keeping only the lin- 
ear terms on 5v n (t) and Sv^(t) we obtain the following 
equation 
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iSv n ,t + kcosq(Sv n+ i + Sv n -i - 2Sv n ) + (A + 2 7 i/ 2 )(<Si>„ + <Sv*) = 0. 



Considering a modulation in the form 

*«»(*)= (^ i(Qn+m) , (23) 

where Q and Q are the wave number and the frequency of 
the linear modulation waves respectively then substitut- 
ing it into l|22|) we obtain the following system of linear 
equations 

ifla + 4fccos<7sin 2 (Q/2)/3 = 0, 
-(4/ccosgsin 2 (Q/2) + 2(A + 2-yv$)v$)a + = 0.(24) 

I 



From the gain equation we can see that the unstaggered 
solution (g = 0) are unstable whenever Az/q + 2 7 j/q > 
and the staggered solutions (q = 7f) are unstable when- 
ever \vq + 2jvq < which leads to the following cases: 

i) q = 0. If A > 0, 7 < then we have MI of the 
unstaggered solution for | 7 | < A/2t/ 2 and if A < 0, 7 > 
then we have MI when | 7 | > \/2v 2 . 

ii) q = 7r. If A > 0, 7 < then we have MI of the 
staggered solution for I7I > A/2i/ 2 and if A < 0, 7 > 
then we have MI when | 7 | < X/2v 2 . 

iii) In contradiction with the continuum case which 
is stable for the set of nonlinearities A < 0, 7 < the 
staggered solution is modulationally instable regardless 
the value of any parameter. However the unstaggered 
solution should satisfy the inequality Ai/ 2 + 2 7 j/q > 2k, 
which is almost satisfied all the time for strongly localized 
modes, in order to have MI. 

In the general case where q can take any value we can 
distinguish two regions of MI as it is shown in Fig. 1 for 
positive k cosq. 

One is fully unstable and the other is conditional means 
it depends on the wave number q of the carrier wave and 
the wave number Q of modulation wave as it is shown in 



(22) 

I 

In order to have a nontrivial solution for the system the 
determinant of the coefficient matrix should be zero that 
leads to the following dispersion law 



(25) 



(26) 

I 

Fig. 2. 

The conditional region represents the effect of the dis- 
creteness on the modulational instability. Clearly we can 
see that if we fix the cubic nonlinearity parameter and 
keep changing the quintic nonlinearity parameter then we 
will have instability, conditional instability and stability 
which shows that the quintic term can lead to collapse of 
MI. The role of quintic term also is very crucial on the 
wave numbers of the carrier and modulation waves lead- 
ing to MI, where increasing 7 in the negative direction 
will shrink the MI region in (q, Q) plane till a critical 
value 7cr . This eliminates all possible chances to get 
MI then MI occurs again in the (q, Q) plane in a mirror 
symmetry with the one occured for 7 < 7 cr . In the con- 
tinuum limit, when Q << 1 and q << 1 equation l|25|l 
reduced to 

n 2 = 2kQ 2 [kQ 2 /2 - (A + 2 1 vl)ul\ (27) 

and coincides with the one obtained in [21J . This MI gain 
does not depend on the wave number q of the carrier wave 
contrary to the discrete case. 

Based on the analytical results, the modulational in- 
stabilities of carrier wave with wave number q modulated 



Q 2 = 8fcsin 2 (Q/2)[2fccos 2 (g) sin 2 (Q/2) - cosg(A + 2 7 ^)^]. 

I 

The plane wave solution (|18|l is modulationally unstable provided that fl 2 is negative and the gain is 
I 

G = 2| sin(Q/2)|W / 2fccosg(A + 2 7 ^)^ 2 - 2fccosgsin 2 (Q/2)). 
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FIG. 1: Modulational instability regions as function of the 
parameters A = j/2kiyo 2 cos q and B = /3/4z/o 2 cos q. with 
cosq > 0, the labelled points are the points discussed in the 
text. 
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FIG. 2: Regions of modulational instability in the plane (q, Q) 
for the conditional stability region in Fig. 1, the labelled 
points are the points discussed in the text. 



10~ 4 . The numerical simulations have been performed 
for a chain of N = 180 units, with periodical boundary 
conditions so that the wave numbers q and Q satisfy the 
relations q — 2'Km/N and Q = 2irM /N, where m and M 
are integers. The other parameters have been chosen to 
be k = 0.1, A = 1 with 7 varying, the amplitude of the 
modulation wave has been chosen small compared to the 
amplitude of the carrier wave as vq = 0.5, Svq — 0.005. 
The stability and instability regions as it was predicted 
by the analytical results has been checked for different 
points , starting with point labelled (1) in Fig.l that cor- 
respond to zero quintic term with q = Q = n/3 and 
q = 7r/3, Q = 27r/3. Our expectation of modulational 
instability has been confirmed as it is shown in Fig. 3 
and independently from q and Q. 



point (3) 

point (4) 




0,2- 



100 200 300 400 500 

time 

FIG. 3: Modulational instability arises at different times at 
point labelled (3) started immediately and at the point la- 
belled (4) the solution is stable for a long time before the 
occurrence of MI. 

The role of discreteness investigated within the condi- 
tional instability region of Fig. 1 with quintic parameter 
7 = —1.2, where we choose two points: the point labelled 
(3) with q — 7r/9, Q — 2ir/3 and the point labelled (4) 
with q = 7r/3, Q = 2tt/3 from Fig. 2 that correspond the 
same point labelled (2) in Fig. 1. The numerical results 
demonstrate the dependence of stability on q and Q as it 
is shown in Fig. 4. 



by small oscillation with wave number Q occurs when the 
right hand side of dispersion equation (|25|l is negative. In 
order to check our analytical results we have performed 
numerical simulations of equation i|14|) using fourth order 
Runge-Kutta scheme with a time step chosen to preserve 
our conserved quantities which are the total energy of the 
system and the number of atoms to accuracy more than 



IV. INTRINSIC LOCALIZED 
MODES(DISCRETE BREATHERS) 

In this section we study the strongly localized solu- 
tions, which are standing and occupying few sites. In the 
case when the coupling between sites are weak the ana- 
lytical approach can be used to find the discrete soliton 
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-point (6) 
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time 



FIG. 4: Modulational instability of the plane wave solution 
occurs and do not occurs at the same point in plane (q, Q) 
with different 7 > y c 7 < y c 
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(breather) solutions. 

We used Page approach [2(| in order to examine the ex- 
istence of the most known localized modes: even, twisted 
and odd modes (see Fig. 5). Twisted and even modes has 
been combined to one mode where both can be described 
by one formula as far as the evolution of one can be other 
after varying the quintic nonlinearity term. 

By considering solution for equation (|14|l as 
v n (t) = y n e( mt ', where v n represents the respec- 
tive amplitudes of a localized mode one can transform 
equation 1|I4[) to a system of algebraic equations depends 
on the topology of the localized mode, hence we derive 
the conditions of existence and then with linearization we 
derive the conditions of stability for each mode as follows. 

Even mode 



FIG. 5: Schematic representation of different kind of odd 
modes (a unstaggered and b staggered), even modes (c un- 
staggered and d staggered)and twisted modes (e unstaggered 
and f staggered). 



Mode 


Existence conditions 


Odd (a) 


A > -7 vl 


Odd (b) 


A < -7 "I 


Even (c) and Twisted (e) 


A > -7 vl 


Even (d) and Twisted (f ) 


A < -7 vl 



TABLE I: Conditions for existences of different kinds of 
modes shown in Fig. 5 and for positive coupling. 



The even/twisted mode defined by its amplitudes 
as v n = Uq (• • • , 0, a, 1, s, sa, 0, • • • ) , where , n = 
±1,±2, as we are concerned with the symmetric 
modes only s = ±1. Varying the signs of s and a 
produces the different topologies of even/twisted modes 
showed in Fig. 5. Taking in consideration the require- 
ments for strongly localized modes which are \a\ <C 1 and 
a n w for n > 2 we derive the dispersion relation and 
the formula of the secondary amplitude respectively as 

u = w E T = k(a + s)+i^(X+jul), (28) 



Hence a set of conditions for existence of strongly lo- 
calized modes can be derived as it is in table I. 



To study the stability of these modes we followed the 
approach developed in 21] by using a linear analysis, 
where we impose a perturbation S n (t) on each non zero 
excitation amplitude, hence the mode's amplitude can be 
written now as 

v n = Vq(- ■ ■ , 0, ra + (5_ 2 , 1 + £-1, s + Si, a + S 2 , 0, ■ ■ ■ ), 

inserting it in l|I4|l and with subsequent linearization we 
got eight-order system of equations. The change of vari- 
ables as Sf = S + j ± 5-j (j=f,2) reduces the system 
to two independent four-order equations systems. Sepa- 
rating the real and imaginary parts of the perturbation 
Sj 1 = 5j r + i$fi and introducing the scaled time r = u) c t, 
we obtain 



7 



dS ± 





(3s — p)ra + 




27^(1 — ra) 
(A + 7^) 



-p)a 





a 


\ 





a 







—a 





1 







-1 





/ 



(30) 



where S ± = {<>i rl S u , S 2r , S 2i ) and p = ±1 stands for the 
symmetric (5^~ and antisymmetric 5~ perturbation respec- 



tively. If we introduce 5^ oc exp(gT) then the eigenvalues 
g of l|3U|) are given by the following equation 



1 



27^g(p -a) 
(A 



+2(p-s) ^1 
2(p-2s) 



•7*8) 
7^o 



2(3 -2ps) c^ + 2(p-s) I 



7^o (g - g) 
A + 7^0 



A + 7^0 
2 7 ^(1 -p + s) 

A + 7^ 2 



2(3 - 2ps) 



271^(1- p + s) 



A + 7^0 



a 4 = 0. 



(31) 



r 



Hence similarly to the case without quintic nonlinear- 
ity |2l| , when the symmetry of the mode coincides with 
the symmetry of perturbation that is p — s the localized 
modes will be stable without any conditions on the non- 



J 



linearity parameters. In the case where the mode and 
the perturbation have different symmetries that p = — s 
then there is a possibility for the mode to be unstable 
with the instability gain 




5s 

T 



r 



7^(1 -2s) 



2(27^ 



7) 



a 



(32) 



In contrast to the case of cubic nonlinearity only, where 
the mode is stable only when as < 0, which means only 
twisted modes ( staggered and unstaggered ) are stable, it 
is clear in the cubic-quintic case that the first coefficient 
of 132(1 leads also to the possibility of stability of even 
modes (staggered and unstaggered) as > 0. However 
taking the existence consideration only the unstaggered 
even mode is stable. The relationships between nonlin- 
earities that control the stability is given in table 2. In 



Mode 


Stability Conditions 


Even (c & d) 


A < -2-yvi 


Twisted (e & f) 


A > -2-yvi 



TABLE II: Conditions for stability of even/twisted modes 
showed in Fig. 1 



the numerical simulations, first we checked the stability 



of the even modes when the symmetry of the perturba- 
tion is the same as the symmetry of the mode. We found 
that the numerical results confirm the analytical predic- 
tions and the modes are stable always. For the case when 
the perturbation and the mode have opposite symmetry 
we checked the validity of analytically derived stability 
regions for the modes. 

The most important result here is that the pre- 
dicted stable even unstaggered mode, which is not 
possible in the case with cubic nonlinearity only, see 
[2l| and references therein, has been demonstrated 
numerically see Fig. 6 for 7 = —0.6 which is less than 
7 = A/2i>o 2 = —0.5. However when 7 = —0.4 which is 
greater than 7 = A/2uo 2 = —0.5, the mode is unstable 
and it is transformed to the odd mode after some time as 
it is shown in Fig. 7. The even twisted modes are shown 
analytically to be stable within two different regions 
and this result has been checked numerically, which is 
in agreement with |2l| for the case of cubic nonlinearity 
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only. 




100 



160 . 



200 



FIG. 6: Evolution of perturbed unstaggered even mode (c) 
A = 1, wo = 1, 7 = —0.6. 




FIG. 7: Evolution of perturbed unstaggered even mode (c) 
A = l,vo = 1,7= -0.4. 

Odd mode 

The odd mode is defined by its amplitude as v n = 
Vq{"- , 0,/3, l,s/3, 0, • • •), where n — 0, ±1,±2, ••• and 
s = ±1 with the requirements for strong localized modes 
are f3 -C 1 and (3 n ss for n > 1 and after performing sim- 
ilar algebraic operations to those done for even/ twisted 
mode we obtain the dispersion relation and the formula 
of the secondary amplitude respectively as . 



P = 



uj = k 2 
k 



1 



"o(* + 7fo) 



2V 



+ A^+ 7 < (33) 
(34) 



which leads to a set of conditions for the existence of 
this mode as it is shown previously in table 1. To study 



the stability of this mode we insert perturbation 8 n (t) on 
each non zero excitation amplitude as 

v n = vq(- ■ ■ , 0, p + 5-!, 1 + 6 , s/3 + 8i, 0, • • • ), 

then substitute it in equation l|14|) and perform lineariza- 
tion procedures, as we did for the even/twisted mode, we 
got a six-order system of equations, again a change of 
variables as = S + i + <5_j reduces the system to four- 
order equations system. Separating the real and imag- 
inary parts of the perturbation 5^ = S +lr + i 5n and 
Sq = 5or + iSoii then introducing the scaled time r = uj c t 
we obtain 



dd od 



t 



\ 





A + 7^n 2 



23 





o p 

-ip 





S od , (35) 



1 J 



where S od = (5$, 6$, S1 d , 5$) . If we introduce 5* cx 
exp(gr) then the eigenvalues g of (|3~5)) are given by the 
following equation 



g i + [l + 4s/?] g 2 + 4 



7^o 2 



A + 7^ 



/3 2 



/3 2 = 0, (36) 



which have four roots and the one that control the sta- 
bility is 



9 = o 



\ 




24P 2 



mxp 2 

\v% + A ' 



(37) 



which shows that the instability gain Re(g) will be a 
nonzero only if 



> P c with p c = ± — 



A + 1V 2 



V8 V A + 3 7 z/2' 

For 7 = 7 C = — A/3^o, P c will be undefined, but we can 
conclude from equation l|36|l that the odd mode will be 
stable whatever the value of p. It is important to men- 
tion here that the value of p c is depend on the interplay 
between the cubic and quintic nonlinearities and satisfies 
the condition for strongly localized modes P <C 1 for a 
range of values of these nonlinearities, in contrast to the 
case of the CDNLSE discussed in [5l| where p c = 1/VS 
which showed that the strongly localized mode is always 
stable because P < p c always. Furthermore for the sta- 
bility of localized odd mode in CQDNLSE the condition 
7 > 7 C = —X/3vq is necessary for focusing cubic non- 
linearity and defocusing quintic nonlinearity. Finally we 
should mention here that for both even/twisted modes 
and odd modes cases, if we substitute 7 = in our re- 
sults, we get back to the cubic case and hence we repro- 
duce the results of [2l| • Numerical results for the case of 
odd modes presented in Figs. 8-9, confirmed the new pre- 
dicted result that is the quintic nonlinearity term can lead 
to an unstable staggered odd mode when —2 < 7 < —1.5, 
this range is chosen to keep the mode strongly localized 
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means the second amplitude 0.1 < (3 < 0.2. However 
the unstaggered odd mode is stable always as it is shown 
in Fig. 10 in agreement with the results for the case of 
cubic nonlinearity only |2l| . 




FIG. 8: Evolution perturbed of staggered odd mode (b) A = 
1, v = 1, 7 = -1.7, p c = 0.140 <j3 = 0.167. 



three-body interactions in the optical lattice. Following 
to the procedure suggested in |jj we derive the lattice 
nonlinear model which describe the GP equation with 
periodic potential. This approximation is reasonable for 
the atomic wave function is approximated by a single 
Wannier function - so called a Wannier soliton. The mod- 
ulational instability in a lattice is the process leading to 
the generation of discrete breathers in the optical lattice. 
We find the regions of instabilities of nonlinear plane 





FIG. 9: Evolution of perturbed staggered odd mode (b) A = 
1, v = 1, 7 = -1.6, (5 C = 0.146 > /3 = 0.143. 

V. CONCLUSION 

In conclusion we have investigated the modulational 
instabilities and discrete breathers in BEC with two- and 



FIG. 10: Evolution of perturbed unstaggered odd mode (a) 
A = l,«o = 1,7 = —0.3. 



wave and conform analytical predictions by the direct 
numerical simulations of the CQDNLS equation. The 
lattice dispersion vary sign in dependence of the quasi- 
momentum value in the Brilloin zone. Thus we can con- 
trol the regions of instability even for the fixed signs of 
two and three body interactions. This circumstance ex- 
tend the possibility for creation of the localized modes 
in BEC in optical potential. We find the expressions for 
the amplitudes and frequencies of different strongly local- 
ized modes in optical lattice such as even, odd, twisted 
stuggered and unstaggered modes. We analyze the stabil- 
ity of these modes and found the stability regions which 
are necessary for the search of these modes in the ex- 
periments with BEC in optical lattice. For the further 
investigations will be interesting to consider the regimes 
beyond of the tight-binding approximation(the deep op- 
tical lattice limit) by accounting the higher order terms 
in the Wannier expansion. In particular it concerns the 
extensions of CQDNLSE involving the discrete nonlinear 
models with nonlocal nonlinearities[7l I22I I23I . 
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